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We predict an unconventional spin-transfer torque (STT) acting on the magnetization of a free 
ferromagnetic (F) layer within N|TI|F vertical heterostructures which originates from strong spin- 
orbit coupling (SOC) on the surface of a three-dimensional topological insulator (TI), as well as from 
charge current becoming spin-polarized in the direction of transport as it flows from the normal metal 
(N) across the bulk of the TI slab. Unlike conventional STT in symmetric F^|I|F magnetic tunnel 
junctions, where only the in-plane STT component is non-zero in the linear response, both the 
in-plane and perpendicular torque are sizable in N|TI|F junctions while not requiring fixed F^ layer 
as spin-polarizer which is advantageous for spintronic applications. Using the nonequilibrium Born- 
Oppenheimer treatment of interaction between fast conduction electrons and slow magnetization, 
we derive a general Keldysh Green function-based STT formula which makes it possible to analyze 
torque in the presence of SOC either in the bulk or at the interface of the free F layer. 



PACS numbers: 72.25.Mk, 75.70.Tj, 85.75.-d, 72.10.Bg 

The spin-transfer torque (STT) is a phenomenon in 
which spin current of large enough density injected into 
a ferromagnetic (F) layer either switches its magnetiza- 
tion from one static configuration to another or gener- 
ates a dynamical situation with steady-state precessing 
magnetization \^ 2^. The origin of STT is absorption 
of itinerant flow of angular momentum components nor- 
mal to the magnetization direction. It represents one 
of the central phenomena of the second-generation spin- 
tronics, focused on manipulation of coherent spin states, 
since reduction of current densities (currently of the or- 
der 10^-10^ A/cm^) required for STT-based magnetiza- 
tion switching is expected to bring commercially viable 
magnetic random access memories [3 . The rich nonequi- 
librium physics [4 arising in the interplay of spin currents 
carried by fast conduction electrons and collective mag- 
netization dynamics, viewed as the slow classical degree 
of freedom, is of great fundamental interest. 

The early phenomenological explanations [5] of STT 
in noncollinear ferromagnetic metal circuits have been 
followed by more microscopic theories j6HTT]. which are 
often combined with first-principles input about real ma- 
terials [BHS] . These theories have been focused on devices 
with no spin-orbit coupling (SOC) where STT is directly 
connected to the divergence of spin current as a conse- 
quence of the conservation of total spin. Thus, STT vec- 
tor can be obtained simply from the local spin current 
at the N|F or I|F interface (N-normal metal, Tinsulating 
barrier) within F'|N|F spin valves or F'|I|F magnetic tun- 
nel junctions (MTJs). Such local spin currents are typ- 
ically computed using the Landauer-Biitikker scattering 
approach [6l |9] or the nonequilibrium Green function 
(NEGF) formalism [IlElIIO]. 

However, this methodology is inapplicable to junctions 




FIG. 1: (Color online) Schematic view of the topological- 
insulator-based vertical heterostructure operated by spin- 
transfer torque. The junction contains a single free F layer of 
finite thickness with magnetization m, which is attached to 
the right semi-infinite N lead. Unpolarized charge current is 
injected from the left semi-infinite N lead across finite thick- 
ness TI slab. We assume that each layer is composed of atomic 
monolayers modeled on an infinite square tight-binding lat- 
tice. 



with strong SOC, which has recently ignited intense the- 
oretical efforts [H [T2HT5] to devise approaches for effi- 
cient computation of STT in the presence of spin non- 
conserving interactions. For example, SOC can be intro- 
duced into the device by the bulk ferromagnets (as in F 
layers based on ferromagnetic semiconductors [2( [T3l IT4] ) , 
or due to the Rashba SOC at the F|I interface in devices 
with structural inversion asymmetry [2]. The impor- 
tance of the latter for potential applications was demon- 
strated in the very recent experiment [16 measuring "SO 
torque" [21 |12] in Pt|Co|A10a, semi-MTJ where charge 
current flows within the plane of a Co layer. 

Concurrently, the recent discovery [17 of three- 
dimensional (3D) topological insulators (TIs), which pos- 
sess a usual band gap in the bulk while hosting metal- 
lic surfaces whose massless Dirac electrons have spins 
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FIG. 2: (Color online) (a) The angular dependence of STT 
components, T|| = r||mx (mxe^) and = r^mxe^, acting 
on the free-layer magnetization m in N|TI|F semi-MTJ shown 
in Fig.[l] (b) The STT components, = r||m x (m x m') 
and — r±m x m^, acting on the free- layer magnetiza- 
tion m = (sin^cos(/),sin^sin(/),cos^) in conventional F^|I|F 
symmetric MTJ where magnetization of the reference layer 
F' is fixed at m' = e^. (c) The STT components in N|I|F 
semi-MTJ, defined in the same fashion as in panel (a), with 
the Rashba SOC of strength an /2a = 0.1 eV located on the 
monolayer of F which is in contact with I barrier, (d) The 
angular dependence of conductances for N|TI|F, F^|I|F and 
N|I|F junctions. The bias voltage Vb in all panels is suffi- 
ciently small to ensure the linear-response regime. 



locked with their momenta due to the strong Rashba-type 
SOC, has led to theoretical proposals to employ such ex- 
otic states of matter for STT applications. For example, 
magnetization of a ferromagnetic film with perpendicular 
anisotropy deposited on TI surface could be switched by 
interfacial quantum Hall current [18 . 

However, very little is known about the device geome- 
tries [19| in which charge and spin currents are perpen- 
dicular to the surface of TI and their potential for ap- 
plications in conventional [4 vertical MTJ setups. In 
general, vertical Tl-based heterostructures would exploit 
strong interfacial SO coupling without requiring [i9\ [20] 
perfectly insulating bulk whose unintentional doping in 
present experiments obscures [2T topological properties 
anticipated for lateral transport along the TI surface. 

In this Letter, we derive an efficient NEGF (in Keldysh 
formulation) -based formula which makes it possible to 
analyze STT in the presence of arbitrary SOC within the 
device. Unlike the recent formulas [131 IE] developed to 
treat SOC eflFects on STT in the linear-response regime, 
ours can handle torque driven by finite bias voltage (re- 
quired to reach sufficient current density in MTJs [4 ), 
and it can also be easily combined with density functional 
theory (DFT) through the NEGF-DFT formalism [8] [22]. 
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FIG. 3: (Color online) The spin-polarization vector 
P°-* = (0, 0, Pr') of current [23] in the right N lead of N|TI|N 
junction as a function of the thickness dxi of the 3D TI layer 
after unpolarized charge current is injected from the left N 
lead. 



This STT formula is then applied to predict unusual 
features, shown in Fig. [2|a), of the in-plane Ty torque 
emerging in Tl-based semi-MTJ illustrated in Fig. [l] in 
the absence of any external spin-polarizer. For conven- 
tional F'|I|F MTJs, where F' reference layer with magne- 
tization plays the role of an external spin-polarizer, 
it is customary to analyze the in-plane (originally con- 
sidered by Slonczewski [5 ) and perpendicular (or out-of- 
plane) torque components, T = Ty + T^. The in-plane 
component Ty = r||m x (m x m') is purely nonequilib- 
rium and competes with the damping. The perpendicular 
torque = r±m x arises from spin reorientation at 
the interfaces and possesses both equilibrium (i.e., inter- 
layer exchange coupling) and nonequilibrium components 
which act like an effective magnetic field on the magneti- 
zation m of the free F layer [11 . While component is 
vanishingly small in metallic spin valves [3 [9] , it can be 
substantial [4 in MTJs due to the momentum filtering 
imposed by the tunnel barrier [IQl [TT]. 

To elucidate the effect of TI slab on unpolarized charge 
current injected from the left N lead we analyze the spin 
density matrix Pgp*^ = |(1 + P°^* • &) for an ensemble of 
outgoing spin in the right N lead of N|TI|N junction. The 
expression for Pgp*^, or equivalent spin-polarization vec- 
tor po^*, was derived as Eq. (10) in Ref. [23] in terms 
of the scattering matrix of the device. Its evaluation for 
N|TI|N junction is plotted in Fig. |3j which shows how 
TI slab polarizes the incoming current in the direction 
of transport with P^^* = (0, 0, 0.5). The polarizing ef- 
fect of the TI slab comes from the effective momentum- 
dependent magnetic field along the z-axis [encoded by 
the Fs term in the TI Hamiltonian in Eq. (|3| discussed 
below], which requires sufficient thickness of the slab to 
manifest itself as well as that the Fermi energy of the 
device is within the bulk gap of the TI slab. The 
spin-polarization of charge current induced by its fiow 
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through a finite-size region with SOC has been discussed 
previously for low- dimensional systems (such as the two- 
dimensional electron gas with the Rashba SOC [24]). 
Due to the constraints imposed by the time-reversal in- 
variance, such SOC-induced polarization cannot [24 be 
detected via current or voltage measurement on stan- 
dard two-terminal ferromagnetic circuits, as exemplified 
by Fig.|2|d) where conductance of N|TI|F junction is the 
same for m || and m ^ configurations. 

Following this analysis, the meaning of torque compo- 
nents plotted in Fig. |2|a) for semi-MTJ is defined by 

T = T|| + Tx = r||m x (m x e^) + T±m x e^. (1) 

In fact, the same definition of torque components is ap- 
plicable [15 to N|I|F semi-MTJ with the strong Rashba 
SOC, aii{& X k||) -e^, at the I|F interface [2, 16 . In that 
case, T|| and components plotted in Fig. [2|c) are 
driven purely by the surface Rashba SOC, which is the 
second order effect oc characterized by torque asym- 
metry [15] around the stable magnetic state 6 = 90°. In 
contrast, Ty in N|TI|F semi-MTJ are non-zero at 6> = 90° 
due to the summation of asymmetric contribution driven 
by the strong SOC on the surface of TI layer and symmet- 
ric [akin to torque in the usual MTJ shown in Fig. [2|b)] 
contribution generated by the conventional STT due to 
spin-polarization of the current passing through the bulk 
of the TI layer. The existence of Ty and for N|I|F or 
N|TI|F semi-MTJs makes this type of STT quite different 
from SO torques [2j [121 [H] that act only as an effective 
magnetic field which can induce switching but not pre- 
cession of the magnetization in the free F layer [15] . 

The exploitation [3 of STT in conventional F'|I|F 
MTJs demands a compromise between large current den- 
sity (requiring low junction resistance to avoid damage) 
and readability (requiring large magnetoresistance). In 
addition, optimization of spin polarization across the 
junction, stabilization of the fixed magnetization of the 
reference F' layer and minimization of stray fields de- 
mands complex stacking structure (involving typically 
more than ten different layers [ln4j). On the other hand, 
our semi-MTJ requires only one F layer. Furthermore, 
Figs.[2];a),(b) show that Ty in N|TI|F semi-MTJ is com- 
parable to the one in F'|I|F MTJ tuned (via the on-site 
potential in the I layer) to have similar resistance. The 
angular dependence of conductances for N|TI|F, N|I|F, 
and F'|I|F junctions are compared in Fig. [2|d). 

Now we turn to the detailed explanation of our formal- 
ism. The semi-MTJ in Fig. [l] is modeled on a cubic lat- 
tice, with lattice constant a and unit area = □, where 
monolayers of different materials (N, F, TI) are infinite 
in the transverse x?/-direction. The TI layer has thick- 
ness (ixi = 5 and the free F layer has thickness dp = 70 
monolayers. The F and N layer are described by a tight- 



binding Hamiltonian with a single 5-orbital per site 

n,(T(T',k|| 



7 Yl (cL,k||Cn+l,a,k|| +H.C.). 
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The operators cj^^ (cna) create (annihilate) electron with 
spin a on monolayer n with transverse momentum ky 
within the monolayer. The in- monolayer kinetic energy 
is ^n,k|| = — 2^ {cos ky a -\- cos kz a) ^ whose effect is equiva- 
lent to an increase in the on-site energy, and the nearest 
neighbor hopping is 7 = 1.0 eV. The coupling of itin- 
erant electrons to collective magnetization dynamics is 
described through the material-dependent exchange po- 
tential = 1.0 eV (A^ = within semi-infinite ideal N 
leads), where & = {d-x^d-y^&z) is the vector of the Pauli 
matrices and [(Jc^Jcrcr' denotes the Pauli matrix elements. 

The minimal model for the slab of 3D TI, such as 
Bi2Se3, is the effective tight-binding Hamiltonian with 
four orbitals per site [26] : 



TI 



2a ^ 



Cn+l,k,| + H.C. 



CI + (i(ky )ro H (Fi sin k^a + T2 sin kya) 



^n,k|| 



It yields the correct gap size in the bulk and surface dis- 
persion while reducing to the continuum k • p Hamilto- 
nian in the small k limit. Here c = (c+^, c+|, c_^, c_|)-^ 
annihilates electron in different orbitals, <i(k)y = M — 
2B/a^ + 25(cos k^a + cos kya - 2)/a^ {i = 0, 1, 2, 3) 
are 4x4 Dirac matrices and 1 is the unit matrix of the 
same size. The numerical values of parameters are chosen 
as: M = 0.3 eV; A = 0.5 aeV; and B = 0.25 a^eV. The 
Fermi energy of the whole device is set at Ejr = 3.1 eV, 
and the bottom of the band of the TI layer is shifted by 
C = 3.0 eV. 

The hopping 7c = 0.25 eV between the sites of metallic 
F or N layer and the TI layer is chosen to ensure that the 
Dirac cone on the surface of TI is not distorted [19t [20] 
by the penetration of evanescent modes from these neigh- 
boring metallic layers. The weak F to TI coupling can be 
achieved by growing an ultrathin layer of a conventional 
band insulator, such as In2Se3 with large bandgap and 
good chemical and structural compatibility with Bi2Se3 
where sharp heterointerfaces have already been demon- 
strated by molecular-beam epitaxy growth [27 . We as- 
sume that such layer of sufficient thickness is present and 
it suppresses the magnetic proximity effect, i.e., A^ = 
on the TI monolayer (denoted as F|TI interface in Fig.jl]) 
closest to the F layer. 

Using the operators cj^^ (c^cr) which create (annihilate) 
electron with spin a on monolayer n, we can introduce the 
two fundamental objects [25] of the NEGF formalism — 
theretarded G^f {t,t') = -ie{t-t'){{cnAt),cl^,{t')}) 
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and the lesser G^;^/ (t,t') = i{ci>^,{t')cna(t)) GF that 
describe the density of available quantum states and how 
electrons occupy those states, respectively. Here (...) 
denotes the nonequilibrium statistical average In 
stationary problems, G^ and G^ depend only on the time 
difference t — t' or energy E after the Fourier transform. 

In the absence of SOC, one can obtain STT in F'|N|F 
spin valves or F'|I|F MTJs by computing [10 the vector 
of spin current between two neighboring monolayers n 
and n + 1 coupled by the hopping parameter 7: 



n,n+l 



Air 



J 



di?dk|| Tr^ [a{G<:,Z 



)]■ (3) 



The integration over ky is required because of the device 
translational invariance in the transverse direction. Since 
for conserved spin current, the monolayer-resolved [9] 
STT is given by = -V-I^ = In-i,n-In,n+i. ^^e total 
torque on the free F layer is [10 , T = Xlv=o(-'-n-i,n ~ 
In,n+i) = 1-1,0 - 1^,00 = 1-1,0- Here the subscripts 
-1 and refer to the last monolayer of the N or I bar- 
rier and the first monolayer of the F layer, respectively. 
In the multilayers with SOC such as those in Fig. [l] this 
straightforward NEGF strategy to get STT becomes inap- 
plicable since spin current will not decay (i.e., I^^^o 7^ 0) 
if SOC is present in the bulk of the free F layer fM]. Also, 
spin current across the interface I^^l q insufficient to get 
STT if strong SOC is present directly at the interface. 

To derive a general NEGF-based expression for a mean 
current-induced force, we start by assuming that the de- 
vice Hamiltonian depends on a variable q which corre- 
sponds to slow collective (i.e., "mechanical") degrees of 
freedom. The expectation value of the corresponding 
canonical force Q = —dH/dq is obtained using the den- 
sity matrix p = J dE G^{E^q): 



+00 

= -— f 

27ri J " 



dETv 



dH 
dq 



dH 
dq 



(4) 



where G^{E^ q) is adiabatic GF obtained for a frozen-in- 
time variable q. By exchanging the derivative between 
the Hamiltonian and G<{E,q), Q = -d{HG<)/dq ^ 
{HdG^ /dq), and by using the standard equations for the 
retarded and lesser GFs [25], G'^iE) = [EI - H - t"^]'^ 
and G<{E) = G'^{E)t<{E)G''{E), we finally obtain 



s<G^r 



dq 



(5) 



The advanced GF is given by G^ = [G^t and / 
is the unit operator. For devices where electron- 
electron or electron-phonon interactions can be ne- 
glected, S^(£;) = T.pt^p{E - eVp) is the sum of re- 
tarded self-energies due to the coupling to semi-infinite 
ideal (F or N) leads p = L,R, TpiE - eVp) = ilt^p^E - 



eVp) — T,p{E — eVp)] is the level broadening operator, 
and T^^{E) = ^^ifp{E)Tp{E — eVp) is the lesser self- 
energy [25). The junction is biased by the voltage eVh = 
eVL — gVr and fp{E) = f{E — eVp) is the Fermi func- 
tion of the macroscopic reservoir to which the lead p is 
assumed to be attached at infinity. 

The expression Eq. ([5| is the central formula of our 
formalism. We note that this STT formula is akin 
to the mean value of time-averaged force in nonequi- 
librium Born-Oppenheimer approach [28 to current- 
induced forces exerted by conduction electrons on ions 
in nano junctions or mechanical degrees of freedom in na- 
noelectromechanical systems whose collective modes are 
slow compared to electronic time scales. The applica- 
tion of Eq. ([5| to get {a = x,y, z) component of the 
STT vector acting on the magnetization of the free F 
layer within N|TI|F junction proceeds by first comput- 
ing G'^{E) for the device described by the Hamiltonian 
H = Hti + Hf- In the second step, the Hamiltonian of 
the F layer is modified 



^1 



Hp + q 4a,k|| [ea • (m X (t)] 



P (6) 



and G^{E)[H^] is computed for the new Hamiltonian 
i^^ = ^Ti + Hf- This allows us to obtain dG'^/dq ^ 
(G^[^^] - G'^[H])/q where we use q = lO'^ as the in- 
finitesimal. The derivative dG'^ /dq plugged into Eq. ^ 
yields Q = Tc,. 

Equation ^ includes both the equilibrium 
T^{Vb = 0) jSl TlQlfn] and experimentally measured [4] 
nonequilibrium T _\_{Vi)) — T _\_(Vi) = 0) contribution to 
T^. The linear-response contribution can be extracted 
by expanding the density matrix p to first order in the 
applied bias voltage Vb and by subtracting the purely 
equilibrium term peq = J dElmGQ{E)f{E): 



dq ^'^'^ '5/^ 




dETi 



dG^dH^ _ dGl dtl 
dq dVp dq dE 



•(7) 



Here Gq{E) is the retarded GF at zero bias voltage 
and we assume zero temperature. The second sum in 
Eq. ^ is non-zero only for cx Vb where the inte- 
gration over the Fermi sea is necessary to ensure the 
gauge invariance (i.e., invariance under a global poten- 
tial shift Vp ^ Vp^ U) of Tx plotted in Fig. [g) N ote 
that (X Vb component is identically zero [4 ^1 0[ [TT] 
in symmetric F'|I|F MTJs, as confirmed by Fig. [2|b) us- 
ing our general Eq. ^ rather than the usual [Uj special 
choice Vl = —Vr = —Vb/2 applicable only to MTJs with 
identical F' and F layers. 

We conclude by noting that one of the key experimen- 
tal issues for STT in conventional MTJs is its control 
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via finite bias voltage [4] . While our Eq. ([5| intrinsically 
takes into account finite bias voltage, the effective Hamil- 
tonian Eq. (|3| is too crude to describe the band structure 
of a real 3D TI material necessary for such calculations. 
Similarly, computation of cx 14 for N|TI|F semi-MTJ 
requires integration in the second term in Eq. ([7| over 
the whole energy band so that our result for is also 
crude. Reliable integration over energy or finite bias cal- 
culations necessitate coupling of Eq. ([5| to NEGF-DFT 
formalism ^2 capture band structure and interface 
reconstruction, as well as self-consistent charge and spin 
densities across the junction. 
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